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Abstract 

It is common to model the random errors in a classical measurement 
by the normal (Gaussian) distribution, because of the central limit theo- 
rem. In the quantum theory, the analogous hypothesis is that the matrix 
elements of the error in an observable are distributed normally. We obtain 
the probability distribution this implies for the outcome of a measurement, 
exactly for the case of traceless 2x2 matrices and in the steepest descent 
approximation in general. Due to the phenomenon of 'level repulsion', the 
probability distributions obtained are quite different from the Gaussian. 
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In classical physics, there is a well-established 'standard model' of random 
errors in measurement, the Gaussian or normal distribution 1 : the error is a 
sum of a large number of more or less independent random contributions, so 
that the central limit theorem assures us that it is a Gaussian. There is not 
as yet a similar theory of errors in quantum measurements. In the quantum 
theory, an observable A is represented by a hermitian matrix. A simple model 
for the error would again be that it is the sum of a large number of independent 
random additions to A. If each matrix element is independent of the others 
(except for the condition of hermiticity), the error would be described by a 
Gaussian random matrix B added to the observable. A measurement of A in 
the presence of this error would yield an eigenvalue of the sum R = A + B 
instead of A. 

For such a random matrix R, the analogue of the eigenvalue problem is to ask 
for the probability distribution of the eigenvalues. Wigner solved such problems 
in the context of nuclear physics EJ , followed by fundamental contributions 
by Dyson [7], Mehta [8], Itzykson and Zuber ^3] an d many others. By now this 
theory has an extensive literature in physics [B] and mathematics E3] • To 
make this paper more accessible, we will solve the problem of determining the 
probability distribution of the eigenvalue for a simple example by elementary 
means first. Then we turn to the more general case, using results from the 
literature on random matrix theory. 

Let us recall in some more detail the theory of errors in classical physics. An 
observable is a function A : M — > R from the phase space M taking real values. 
The state of a classical system is given by a point £ G M in the phase space. 
Ideally, a measurement of the observable A on a system in this state will yield 
the value A(£).The standard model of errors in classical measurement is that 
there are a large number of small, more or less independent, random corrections 
Bx, ■ ■ ■ Bm which add to this observable: 

M 

R = A + Y,B r (1) 

k=l 

Irrespective (upto mild assumptions) of the distribution of Bx ■ ■ ■ Bm, the sum 
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will, in the limit M — > oo will tend to a Gaussian distribution ( the central limit 
theorem) JJ]: 

P{B) = e- cB \ (2) 

We can assume that the mean of this distribution is zero since otherwise it 
can be absorbed into the definition of A. ( In any case such systematic errors 
cannot be analyzed by statistical methods.) The variance is usually assumed to 
be the same at all points in the phase-space so that c is a constant. Thus we 
model the outcome of the measurement by a Gaussian random variable whose 
mean is the 'true' value a — A(£) of the observable and the standard deviation 
is a measure of the size of the error: 

PA^ae^H 1 . ( 3 ) 

There is by now a well-developed sampling theory on how best to estimate this 
mean a and standard deviation a from repeated measurements of R ■ 

If the state of the system is not known exactly, there is a probability den- 
sity function on the phase space p : M — > R which determines the 'instrinsic' 
probability distribution of A in the absence of errors: 

P A (a) = J 5(a-A(0)p(m (4) 

If we add in the error, we get the convolution of this with the Gaussian: 

p A (x) oc J P A {x - B)P(B)dB (5) 

Thus there are two separate ways in which statistical considerations enter 
classical error analysis: the classical observable can have an intrinsic randomness 
because the state of the system is not completely known; and there can be a 
random error added to the observable. In the quantum theory, in addition there 
is another, more fundamental source of randomness: the uncertainty principle. 
Thus a proper theory of errors has to take into account all three sources of 
randomness. 

In the quantum theory, an observable [31 E] is represented by a self-adjoint 
operator on the Hilbert space of states. Let us consider a hermitean matrix A of 
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finite dimension N with eigenvalues a\ , ■ ■ ■ ajv and corresponding eigenstates 
\u\ > ■ ■ ■ \un >• Ideally, if a measurement of the observable represented by A 
is made on a system known to be in state \ip >, the outcome will be one of 
the eigenvalues at, with probability | < Uk\ip > | 2 - Another way to state this 
result (convenient in the following) is that the probability density function of 
this random variable is 

P A (x) =< i/>\5(x - A)\1> > . (6) 

More generally if the state of the system is only partially known, we have 
a mixed state described by a density matrix p which is a positive hermitean 
matrix with tr p = 1. We have then the 'intrinsic' probability distribution 

P A (x) = tvp5(x-A). (7) 

This formula takes into account two of the sources of randomnes mentioned 
above: the fundamental uncertainty of quantum mechanics as well as our pos- 
sible lack of knowledge of the state of the system. 

However, this is still an ideal situation that is never realized in practice: there 
is always some error in the measurement which is a third source of randomness. 
It is useful to have theory of errors in quantum measurement analogous to 
the classical theory above, particularly in view of current interest in quantum 
information theory. 

As an elementary example, suppose the observable is a traceless 2x2 matrix 
„ fa \ 

A = . We can think of it as the energy of a spin half particle in a 

VO -a) 

constant external magnetic field along the z-axis. There are different sorts of 
errors that can affect this energy. The simplest possibility is that the magnetic 
field fluctuates due, for example, to thermal effects in the currents producing 
it. We can model this as a sum of small independent additive corrections to the 
magnetic field. By the central limit theorem, this additional magnetic field B 
can be represented as a Gaussian random variable. 

R = A + a B. (8) 
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The mean of B can be assumed to be zero. Otherwise it can be absorbed into 
the definition of the fixed magnetic field; after a rotation we can reduce it to the 
above diagonal form anyway. It is reasonable, and simplest, to postulate that 
the error is rotation invariant;i.e., all the components have the same variance: 

P(B)oce- cB \ (9) 

We can also think of this as the Boltzmann distribution for the fluctuating 
magnetic field, as the energy of a magnetic field is proportional to B 2 . This 
means that the matrix elements of R are also Gaussian random variables, with 
mean A. 

The eigenvalues of R = A + er ■ B will be ±r, where r is a real random 
variable. Calculating its distribution from that of R, 

P A {R)ote-i ctI ^- A f, (10) 

is now the analogue of solving the eigenvalue problem. We will first consider the 
case where all the states are equally probable: the density matrix is a multiple 
of the identity. In effect we have to average over all fluctuations that change 
R without affecting its eigenvalues. Representing R = a ■ R as a vector and 
transforming to spherical polar co-ordinates, we can perform the average over 
random direction of the vector R: 

p A (r) cx J tr l -5{r-R)P A (R)dR 

cx r 2 J e- c( - r2+2arcos8 hmede (xrsmh[2acr}e- cr2 . (11) 

Recalling that the observed eigenvalue can take also negative values, we nor- 
malize this distribution to get 

p A (r) = ^^rsmh[2acr]e- c{r2+a2 \ (12) 

We can also write this in a way that explicity displays the peaks at r = ±a: 

p A (r) = ^\j^r [e- c ( r - a ) 2 - e^^) 2 ] (13) 

In Fig. 1 we plot this probability distribution. 
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In summary, the quantum mechanical error is modelled by a Gaussian for 
the matrix elements which leads to a markedly distribution from the Gaussian 
for the observed eigenvalue. There is a second order zero for the distribution at 
the origin, a consequence of the 'level repulsion' of the eigenvalues of a random 
matrix. Note that the peaks are displaced outwards from the 'true' eigenvalues 
±1 due to this level repulsion. (In the figures we have assumed an unrealistically 
large error to illustrate the phenomena better.) 

The phenomenon of level repulsion is well-known in the theory of random 
matrices. The set of all traceless hermitean 2x2 matrices with a given spectrum 
±r is a sphere of radius r. The volume of this sphere shrinks as r — > 0, so that 
it is unlikely that the eigenvalues of a random matrix are close together. We 
can also think of the logarithm of the volume of the set of all matrices with a 
given spectrum as an entropy |11|. The probability distribution above can be 
thought of as a compromise between maximizing this entropy and minimizing 
the energy of a fluctuation in B. 

So far we dealt with the case when all the states of the system are equally 
likely;i.e, when the density matrix is proportional to the identity. The proba- 
bility distribution in general will be given by averaging over the error as well as 
the states weighted by the density matrix: 



In the special case of a two dimensional Hilbert space, we can expand 



The inequality ensures that the density matrix is positive; the trace of p is 
normalized to one. We can again evaluate the integral by passing to spherical 
polar co-ordinates, using the identity 




(14) 



p= - +(T . P: \ p \ < _ 



(15) 




a R 



[5{x — r) — 5(x + r)] 



(16) 



2r 



to get 
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where r) = tr pA is the average of A over states. In Fig. 2 we plot this 
probability density in the extreme case of a system in an eigenstate of A. Note 
that there is a small peak at the 'wrong' eigenvalue, caused by the errors. 

We are ready now to take up the more general case of an N x N hermitean 
matrix A with an additive random error B. We can assume that the mean of B 
is the zero matrix: otherwise, we can redefine A by absorbing this mean matrix 
into it. For simplicity we assume that the matrix elements of B are statistically 
independent of each other ( except for the constraint of hermiticity) . In order for 
this condition to be true in all choices of basis, the joint probablity distribution 
should be a Gaussian jHj: 

P(B) oc£r Cltri ^ C2 ( tr ^ 2 , (18) 

for some positive constants c±, ci. Another justification of the normal distribu- 
tion would come from the central limit theorem: each matrix element of B is 
the superposition of a large number of small errors. 

In the presence of the error, we are really measuring the observable R = 
A + B. The joint probability density function of its matrix elements is 

P A (R) oc e~ Cl tr iA-Af-c 2 [tr R— tr A] 2 (1Q) 

We can write R = UrW where r is a diagonal matrix and U a unitary matrix. 
By averaging over U we can get the joint probability density of the eigenvalues 
of R. In this process we must remember that the Jacobian for transforming 
from the matrix elements of R to r, U is [S] 

dRocA(rfd N rdU, A(r) = *[[(r k - n). (20) 

k<l 

Thus the joint distribution for the eigenvalues of R is 

p A {r u ■ ■ -r N ) cx A(r) 2 e" Cl E k H+ a lh c 4E k r *-£ fc / dUe 2ci tr r(jiM 

(21) 
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The last integral was evaluated by Harish-Chandra ( an d rediscovered by 
Itzykson and Zuber 13 in a more physical context) 

U(N) A(r)A(a) 



Thus 



det e 2cirkai e" Cl E* [ r 2+°2]- c *E fc ^"E^ a <=] 2 



Pi(ri,---r n ) oc =^dete" Cirfcai e" Cl ^a'^ a fe J- c ^Z. fe ^-Z. fc a '=J (23) 
A(a) 

We can also write this result as 

p*(ri,--Tjv) oc ^fl Y sgn ( CT )e- Cl E,h-^J 2 -«[E fc ^-E fc a *] 2 (24) 
Ala) ^— ' 

where the sum is over all permutations of the indicesjl, 2, • • ■ N}. 

fa \ 

For the simplest example -A = we considered earlier, we can write 

\0 -a) 

the above as 

p A (n,r 2 ) oc e -[ c 2+| C i]ki+r-2] 2 ( ri _ r2 ) s inh[2aci(ri - r 2 )]e-^ Cl{ri - r2) \ (25) 

This shows that the sum and difference of eigenvalues are independent random 
variables. The sum is a Gaussian random variable with mean zero. The dif- 
ference has the distribution we derived earlier. Thus c 2 describes a 'classical' 
source of error which shifts both eigenvalues the same way (and is a Gaussian) 
while ci is a 'quantum' error that affects their difference. 

By integrating over all but one eigenvalue we can get the probability distribu- 
tion for the outcome of a single measurement of A. This leads to a complicated 
expression in general but in the steepest descent approximation we can get a 
simpler form: 

N . , 

p A (x) oc^e"^--) 2 J] (26) 

k—1 k^ra 

where c = c\ + c 2 . This probability distribution is peaked near (but not at) 
the eigenvalues a\, • • • a at as one might expect. However even when the density 
matrix is proportional to the identity, not all these peaks are of the same height: 
due to the repulsion of eigenvalues of a random matrix, the extreme values are 
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more probable. Moreover, the peaks of the probability are displaced towards 
the edges due to this repulsion. We plot an example in Fig. 3. 

The spacing between eigenvalues are of order fi. In the limit as Ti — > and 
N — > oo we get back the classical description. It is possible to develop a semi- 
classical theory along the lines of Ref . ^3] . We hope to return to these issues in 
a later publication. 

It would be interesting to verify our predictions experimentally. Perhaps 
Superconducting Quantum Interference Devices ( SQUIDs) or Spin Resonance 
will provide such tests. 

Acknowledgement I thank A. Agarwal,G. Agrawal, J. Eberly, G. Krish- 
naswami, E. Langmann, A.C. Melissinos and C. Stroud for discussions. 
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Figure Captions 



1. Probablity density of the outcome of measuring an observable with A 
eigenvalues ±1 in a mixed state with equal weight for the eigenvalues and 
error parameter c = 3. The peaks are shifted to ±1.15 from ±1. 

2. Probablity density of the same observable in an eigenstate with eigenvalue 
1.0 and error parameter c = 3. The errors cause a small peak at the 'wrong' 
eigenvalue. 

3. Probability density of an observable with eigenvalues — 1, 0.5, 1.0 with error 
parameter c = 2. The peaks are shifted to —1.25, 0.36, 1.41 respectively. 
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Figure 2 of S. G. Rajeev 



13 



prob density 




outcome 



-1 .5 1 
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